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We present and compare the predictions of various cosmic-ray Monte Carlo models for the 
energy (dE/drj) and particle (dN/drj) flows in p-p, p-Pb and Pb-Pb collisions at y?s NN = 14, 
8.8, and 5.5 TeV respectively, in the range covered by forward LHC detectors like CASTOR 
or TOTEM (5.2< \rj\ <6.6) and ZDC or LHCf >8.1 for neutrals). 



I. INTRODUCTION 



The origin and nature of cosmic rays (CRs) with energies between 10 eV and the Greisen- 
Zatsepin-Kuzmin (GZK) cutoff at about 10 20 eV, recently measured by the HiRes [l| and Auger 0] 
experiments, remains a central open question in high-energy astrophysics. One key to solving this 
question is the determination of the elemental composition of cosmic rays in this energy range. 
The candidate particles, ranging from protons to nuclei as massive as iron, generate "extended 
air-showers" (EAS) in interactions with air nuclei when entering the Earth's atmosphere. Due to 
their low observed flux (Fig.[TJ left), only indirect (yet complementary) measurements are possible 
using the atmosphere as "calorimeter" . The first method relies on measuring the fluorescence light 
emitted by air molecules excited by the cascade of secondaries. The second one relies on the use of 
either scintillators or water Cerenkov tanks to sample the shower at ground. Recent Auger results 
are consistent with the hypothesis that the highest energy CRs are protons (source correlated with 
Active Galactic Nuclei) [3[. However, fluorescence-based measurements of the shower maximum 
as a function of primary energy favour a mixed composition (Fig. [TJ right). 




FIG. 1: Left: Measured cosmic ray energy spectrum compared to the center-of-mass energy of various 
hadron colliders. Right: Measurements of CRs shower maximum in the atmosphere as a function of 
primary energy for several experiments along with MC predictions for proton and iron primaries. 
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The determination of the primary energy (from the surface detectors alone) and mass (from either 
method) relies on hadronic Monte Carlo (MC) codes which describe the interactions of the primary 
cosmic-ray in the upper atmosphere. The bulk of the primary particle production is dominated 
by forward and soft QCD interactions, modeled commonly in Regge-Gribov-based Q approaches 
with parameters constrained by the existing collider data (Ei a b % 10 15 eV). When extrapolated to 
energies around the GZK-cutoff, the current MCs predict energy and multiplicity flows differing 
by factors as large as three, with significant inconsistencies in the forward region. The coming 
energy frontier for hadron collisions will be reached by the Large Hadron Collider (LHC). The 
LHC will boast a full compliment of detectors in almost the full range of pseudorapidity (Fig. [5]). 
Measurement of forward particle production in p-p, p-Pb, and Pb-Pb collisions at LHC energies 
(Eiab ~ 10 17 eV) will thus provide strong constraints on these models and allow for more reliable 
determinations of the CR energy and composition at the highest energies. 



II. FORWARD DETECTORS AT THE LHC 



All LHC experiments feature a detection coverage at all rapidities without parallel compared 
to previous colliders (Fig. [5]) . ATLAS @ and CMS @, § not only cover the largest pr-f] ranges 
at mid-rapidity for hadrons, e^, 7 and /i's, but they also feature extended instrumentation at 
distances far away from the interaction point (IP). Forward calorimetry is available at ±11 m (the 
FCal and HF hadronic calorimeters), at ±14 m (CMS CASTOR sampling calorimeter) Q, and 
at ±140 m (the Zero-Degree-Calorimeters, ZDCs) [13, Ell- In addition, ATLAS has Roman Pots 
(RPs) at ±220,240 m, and there are advanced plans to install a ne w p roton-tagger system at 420 m 
(FP420) from both the ATLAS and CMS IPs QjJ. Both ALICE [l]| and LHCb [14] have forward 
muon spectrometers in regions, 2 < r\ < 5, not covered by ATLAS or CMS. The TOTEM 
experiment [ll|, sharing IP5 with CMS, features two types of trackers (Tl and T2 telescopes) 
covering 3.1 < \t]\ < 4.7 and 5.2 < \rj\ < 6.5 respectively, pl us p roton-taggers (Roman Pots) at 
±147 and ±220 m. Last but not least, the LHCf experiment [Ty] has installed scintillator/silicon 
calorimeters in the same region of the ATLAS ZDCs, ±140 m away from IP1. 
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FIG. 2: Approximate pt-t] coverage of current (and proposed) detectors at the LHC [f|. 



III. MODEL PREDICTIONS 



The bulk of particle production in high-energy hadronic collisions can still not be calculated 
within first-principles QCD. Phcnomcnological models based on general principles such as 
unitarity and analyticity are often combined with perturbative QCD predictions for high-p^ 
processes to obtain an almost complete description of the final states. Such "QCD-inspired" 
models - like DPMJET 0, EPOS neXus 0, QGSJET 01 and II HO, S3, and SIBYLL [H 
frequently used for simulating cosmic ray cascades - use an unitarization procedure to reconstruct 
amplitudes for exclusive processes and to determine the total and elastic cross sections. Soft 
processes are described within Gribov's Reggeon theory Q and hadrons are produced mainly 
in the fragmentation of color strings. Whereas this Regge-Gribov approximation is applied to 
hadrons as interacting objects in the case of DPMJET, QGSJET and SIBYLL, it is extended to 
include partonic constituents in EPOS and neXus. Other important differences among models 
are the approximations of high parton density effects at small values of x = p parfc ,„/Pp r£ , f£ ,„ and the 
treatment of the hadronic remnants in collisions. SIBYLL and DPMJET use an energy-dependent 
transverse momentum cutoff for minijet production based on the geometric criterion that there 
cannot be more gluons in a hadron than would fit in the given transverse space [23]. QGSJET II 
implements enhanced Pomeron-Pomeron interactions which are equivalent to impact-parameter 
and density-dependent parton saturation for soft processes [24j, achieving a good description of 
RHIC Au-Au data. EPOS, on the other hand, uses the wealth of RHIC data to parametrize the 
low- a; behaviour of the parton densities and the string fragmentation parameters [25| . reaching a 
good reproduction of a wider range of experimental results. 

Figures [31 2] and [5] show for p-p, p-Pb and Pb-Pb respectively, the inclusive multiplicity and 
energy flows predicted by the models (PYTHIA 26] p-p predictions, with tunings as used in [27l ]. 
are also included) for all pseudo-rapidities, as well as the energy deposit in the acceptances 
covered by the CASTOR/T2 (5.2< |r?| <6.6) and ZDC/LHCf (|r?| >8.1 for neutrals) detectors. 
In some cases, differences as large as 60% are observed. The excellent and complementary 
detection capabilities of the various experiments at the LHC will allow for a (fast) validation and 
tuning of the MC hadronic models - already with a first p-p run with very low integrated luminosity. 




FIG. 3: MC predictions for p-p collisions at 14 TeV: energy flow (dE/dr/, left) and energy deposited in the 
CASTOR/T2 acceptance {dE/dr)\ M=B . 2 - 6 .6, right). 




FIG. 4: MC predictions for p-Pb collisions at 8.8 TeV: Inclusive particle multiplicity {dN/drj) and energy 
(dE/drf) densities (top), neutral energy flow at rapidities covered by the ZDC/LHCf detectors (bottom). 
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FIG. 5: MC predictions for Pb-Pb collisions at 5.5 TeV: Inclusive particle multiplicity (dN/drj) and energy 
(dE/drf) densities (top), and xf distributions for 7r°s and neutrons (bottom). 
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